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IHTHODDOTIOS 

The method of characteristics for three dimensional axially 
symmetric bodies was used to determine the velocity distribution 
about the nose of the Corporal E rocket, a rocket projectile, for 
Mach numbers 2 and 4. Prom the velocities the pressure distrib- 
ution was determined and a drag coefficient computed. 

For a starting point the no3e of the projectile was approx- 
imated for a short distance by a cone and the Taylor-fiaccoll’* 
solution to this problem %*as used. This solution gave the angle 
of shockwave and the body. The Sauer graphical-numerical itera- 
tion method was used for the remainder of the solution. 

Preliminary calculations and the work for Mach number S were 
carried out by Dr. H. K. Forster to whom v;e aro indebted for 
instruction and aid with this work. 



•Taylor, 0. I. and Maccoll, J. V. , Proc. Roy. 3oc., A, Vol. 139, 
po 278-311 1933. 
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CHAHACTKRISTIC3* 

The following discussion is for two dimensional flow and is 
taken from w Introduction to Aerodynamics of a Compressible Fluid", 
by H. ft. Liepmann and Allen 2. Puckett. 

I 

The term "characteristic" is related by definition to proper- 
ties of certain partial differential equations. For any second 
order, partial differential equation of the type 



(1) A 



€1 



■h ZB 



<)x hy 



t c 



ILL 

Zy* 



= d-j| 4-rii 

4 * 



two families of curves nay be defined by the total differential 
equations 



(s) Ai. = 6 1. y e* - ac 

A 

where A, B, and C are functions of x and y. These curves (2) are 
known as the characteristics of equation (1). 

These characteristics have certain properties which make them 
useful in an approximate solution process. These properties will 
be listed without proof: 

1) They are invariant under transformation. 

2) A solution Z(x,y) to equation (1) nay be composed of 
separate integral surfaces joining continuously with continuous 
derivatives along a line. This line is a characteristic. This 
oroperty of the branching of integral surfaces along a character- 
istic permits us‘ to regard the line as a wave front along which a 
disturbance is propagated through a field. 



•Liennann, H. V. and Puckett, Allen 2. Introduction to Aerodynamics 
of a Compressible Fluid. t>p 229-232 1947 
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3) Boundary conditions along any line, C, in the field will 
define a solution within a triangle bounded by C and the charac- 
teristics through the end points of C. Boundary conditions along 
a characteristic, however, will not define a solution in any region; 
boundary conditions must bo given along two intersecting character- 
istics of different families. The exact equation of the velocity 
potential is: 



(3) 



2 



M. AT 



h*0 
<SX Jy 



-h (a . 2 ~v l 




vrtiere -U ~ jjjr 'V = ^ 

If a new dependent variable is defined 

F (^,^) - stt a t- <vy - / 



then when a transformation is made to the hodograph plane, an ele- 
ment of the F surface in the (u,v) plane is defined by an element 
of the (x,y) surface. (Legendre Transformation). The transform- 
ation of equation (3) becomes 



(4) 



(* l - v~) (^--V’jr vv 




The characteristics of the equation are 

(5) .d V A4, / 1< r - ’V 1 1 ' )( ^ - 'V ) 

d ( (X x - sy- L J 

This equation can be integrated immediately, w a w having been expressed 
as a function of u and v. Thus in the two dimensional case a charac- 
teristic chart may be constructed and is known as the "Prandtl- 
Busemann Characteristic Diagram. M 

When an attempt is made to use this method in three dimensional 
problems it is found that it cannot be done. A Legendre Transformation 
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leads to a non-linear equation In this case. In general, in three 
dimensional flow, it is not nossihle to construct a characteristic 
diagram. Instead, a separate construction must bo carried out for 
each flow. 
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SAUKR METHOD 



The 3auor solution of supersonic axially symmetric flo is a 
graphical-numerical iteration method. 

7or a non-viscous, steady, isentropic flow of a perfect fluid 
without body forcds, the equation of notion may be obtained by 
eliminating the pressure and density from the force equations and 
the equation of continuity. Transformed to cylindrical coordinates, 
and assuming axial symmetry, this becomes 



<«> o- +0- 



a 



a 



1 



he \+ - o 

/ -a 



Assuming irrotational flo*/ the velocity components may be replaced 
by the proper derivatives of the velocity potential <$> . 

In the axially symmetric case, given a velocity vector, V, one 

/ >rthogonal 

axes in which X 1 is along 

the velocity vector. Then u 1 = V and v* * 0 where u 1 is the velocity 
in the X’ direction and v 1 is the velocity in the I* direction. 

If at a particular point the X 1 , Y’ coordinates are introduced 
equation (6) becomes 



(I - ^ ) 
v aw 



4*' 



4X 



t 



j r 

•A 



- o 



and introducing 0 , for irrotational flow 




The Mach angle is defined as c\ - X where a is the local 

soeed of sound and V the local fluid velocity. A Mach line is a 
line inclined to the direction of fluid flow by an angle equal to 
the ilach angle. 



At a point, where the velocity i3 known two Mach lines nay 
be drawn. In this case where the velocity is along the X' axis, 
one f*ech line will be at an angle and one at an angle froa 
the X* a.xis. (See Pig. 1 a). 




These Mach lines, which may be identified as the character- 
istics of equation (6) define a new set of coordinate axes, and 
From the definition of the Mach angle, sin it is seen that 

substituting this in equation (7) we get 




Fig. (1) b 
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substituting these relations in (8) 

_ .XX g JJ* 
-a. 

Shis equation can be written 



(9) 






X- z 7 1^. ) - A3-~ - * iJ! - 

^ 1 J*? ; d] “T - ' 

X /<)j2 ) - ,^£- - -±rr^u-X , a r* 
f? ^ dr ? ' c/^ 



*? v X 

where p and q indicate the projections of the velocity vector along 

the | and Q directions respectively. ^ 

rienlacing the derivatives of (9) by a quotient of increments 

and clearing fractions the equations become 

J 'f . ? v J? 

0 U o \ 

Jp - J£ <X dr} 

and being the increments in the directions of the Mach lines 
and dp and dq the projections of the velocity increment unon the 



(10) 



Mach lines. 
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PROCEDURE 

The properties of the characteristics enable us to make an 
approximate solution. The mechanical process i3 a simultaneous 
construction of velocities in the hodograph plane and of <lach 
lines in the physical plane. Step by step small areas in the 
physical plane correspond to points in the hodograph plane since 
the flow field in the physical plane is divided into a finite 
number of small regions bounded by Mach lines over which averages 
are taken for conmuting velocities in the hodograph plane. Hence 
the accuracy of the computation will defend on the size of the 
regions bounded by Mach lines. The boundary conditions in the 
physical plane determine the direction of tho velocity in the 
hodograph plane; at the shock wave the magnitude is al®o deter- 
mined. The points in the hodograph plane when found by an iter- 
ation process based on estimates will in turn determine the Mach 
lines in the physical plane. 

The two boundaries which define the boundary conditions are 
the shock wave attached to the nose and the surface of the body. 
Between these boundaries a net of Mach lines is constructed by 
proceeding from point to point until, in the entire area about 
the body and at points along the surface, the velocities are 
determined. This Mach net is shown in Fig. (5), Fig. (6) or Fig, 
(7). 

The angle of the shock wave at the nose is obtained from 
Taylor-Maccoll solution of flow about a cone with the same nose 
angle. For this purpose, defending upon the geometry of the body 
the assumption that the nose is conical will be good at least a 
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small distance back from the vertex. Taylor and Maccoll showed 
that velocity magnitudes and directions are constant in super- 
sonic flow about a cone along radial lines through the vertex. 

This solution furnishes as many joints as needed to start the 
characteristic method solution. These points can be taken close 
enough to the vertex to make the conical flow assumption valid. 

V.'ith the boundaries established and some points in the phys- 
ical plane for which the velocities are known one nay proceed by 
using the Sauer graphical method for computing supersonic axially 
symmetric flow. There are three different geometrical cases for 
which slightly different mrocedures are necessary. The cases are: 

(a) from the known solution in two points in the free stream 
to calculate the solution in a third point in free stream; 

(b) from the known solution in a point in the free stream to 
calculate the solution in a point on the body; 

(c) from the known solution of a point on the shock wave end 
a point in the free stream to calculate the solution in a third 
point on the shock wave. 

‘first, take the case of all points in the free stream, (see 

fig. 2 ). Suppose points P^ and P^ are known: to find t > r _ % is 

estimated in the hodograph plane. Then with mean values, magnitude 

and direction, which are estimated, draw an average J and/p 

line. Their intersection gives F in the physical plane from which 

values are obtained for d J and A. . With these values one can 

calculate dq and dp from 

dej - xulv x d A j 

d ~ ^ ^ d ^ 



Ilodograph Plano Phyaical Plans 
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before. The velocities from successive points along the body fall, 
however, in a smooth curve in the hodograph plane making estimation 
easy with the result that the process generally does not need to 

be repeated. 

Now in obtaining a point on the shock wave from a preceding 
point on the shock wave and a point in the free stream it must be 
noted that the angle of the shock wave is not known. The shock 
wave is conical only for flow about a cone. For a curved profile 
the shock wave is also a curve which for this method is approxi- 
mated by a series of straight lines, the angles of which must be 
found. To do this the Busemann shock polar is used in the hodo- 
graph plane. The shock polar graphically expresses the relation 
between the velocity direction and magnitude of the downstream 
flow and the angle of the shock wave for. a given free stream vel- 
ocity. Referring to Fig. (4) suppose P^ and are known; is 
to be found. Since P 3 is on the shock wave it is also on the 3 hock 
polar along with PJ which is known. J59tinate Fj on *ke shock polar. 
The angle of the shock wave in the physical -Diane is the complement 
of the mean of Z OAP^ and ^ OAP^. lino P^F-r is obtained as 

in the free stream case. The intersection of this ^ line with the 
shock wave gives the values necessary to calculate do. The con- 
struction in the hodograph plane is made giving the -point at the 
intersection of the shock polar with the perpendicular to the J line 
from P^. Since the shock wave and the line from are nearly 
parallel it is often necessary to repeat the procedure. 
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Then stepping off dq in the 1^ 13 direction and dp in the J ^ 
direction perpendiculars are erected as shown in Fig. (2). These 
intersect in a point which will he very close or right on the 
estimated F^. If this is not within the desired accuracy the 
process is repeated until the iteration 3 hovs no change. Table I 
shows several examples of calculations made in obtaining the values 
necessary to construct Fig. (2). 

For the case of obtaining a point on the surface of the body 
from a ooint in the free stream the procedure differs slightly. 
Referring to Fig. (3) suppose P^ is known; ?2 found. A 

J line, whose direction is estimated, is drawn from until it 
intersects the boundary at B. Since the flow is tangent to the 
boundary at all points on the boundary the direction of the flow 
velocity at point B is known, and an estimate of the magnitude of 
the flow velocity is made. From the mean between this estimated 
velocity and the velocity at a Mach angle is calculated. Sim- 
ilarly a mean value of the angle of flow direction (between P^ 
and point B on the surface) is calculated. The difference between 
these two angles gives a new J line intersecting the boundary at 
jF'g. The points F^ and ^Pg give the values necessary to calculate 
dq as in the free stream case. As before the construction is made 
in the hodograph plane and the perpendicular at a distance dq from 
along the line intersects *2 l p p» From the mean values of 
and draw the corresponding J line from P^ in the physical 
plane. This should intersect the boundary at not, one 

takes the tangent to the boundary at the new point and proceeds as 



Holograph Plano Physical Plano 
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Hodograph Plana Physical Pl^ne 
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So ?8 
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Then* 



Cp = * 



7 % 



Z 1# a function of the strength of the shock wave and thus of the 
shock wave angle, which is known from the Taylor-Maccoll solution 
(the angle of the shock wave is taken at the nose), and of the 
free stream Mach number. 

Given and c< (see figure below) i? is computed as follows*? 



2 . 




c<r° 



C (r-0 +(n-i)(*)J(r-i) 



i r 

( r- 0 i ( r-t-i ) (jj)_ 

(<r\) f 



3. 



•Durand* W. P. , Aerodynamic Theory, Volume III, Durand Reprint, 
1943, pp 238-240. 
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CALCULATION of C and c 

P D 



After the calculation of the flow velocities is conpletod, 
all of the points of the Mach net "between the shock wave and the 
surface of the body aro tabulated, l.e., the geometrical position 
of the point and the stream velocity at the point, (see Table II ), 
In the calculation of the drag of the nose section, however, only 
the points on the surface of the missile are needed. Given a num- 
ber of points on the surface and the velocities at these points 
the pressure coefficient, C^, is computed for each of these points 
and a plot is made of G p vs. X. (see Fig. (8)). From this curve 
values of 0 ^ are obtained for regular increments of X and then by 
carrying out the appropriate integration (in this case a graphical 
integration was used; see Fig. (10)) the drag coefficient, C^, can 
be obtained for the nose of the missile. 

The equations used in the actual computations are discussed 

» 

below! 



In general 



c, s ^ = t - : 



4 .ovv 1 
2 - / 



and 



fO tv 



- \_P w 1 _ Yjh 

2 * */> ~t~ < 11 1 



where Is the free stream Mach Number, 



c 



r 



£ - £ 



A 



* 



'XJz— « - 



So 



M, 



Across the shock wave, however, there is a drop in stagnation 



pressure. This drop is defined by a = V^/v 0 where P3 corresponds 
to stagnation conditions behind the shock wave and p Q corresponds 
to free stream stagnation conditions. 
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4. 




/ ^ 



Y-l 



y k & 



[( r -0 + (r+ ')(■%)] 








Hi, the free atrean Mach Humber is known from the Taylor-Maccoll 



solution of flow around a cone. H* = 



local velocity 



3 critical speed of sound a 
is known at various points on the surface of the body and 



w 

x- 



( rt s ~) 



or 




n. 



.2 t1< 



(rtt)-(r-t) Mj* 



where M g i9 the Mach number at any point on the surface of the 
missile. Substituting for M in the equation for p./pJ 



or 



P- a 
r 



t = 

h 



r 



I - 



x=l ( 




* c 


(■ 


(rti)-'(r-i) h 


(Jzll 


l 

*s J 


X 

r-/ 


(n<) 





-i X. 



Substituting these quantities into the expression for C : 

P- - fu p 



C 



P 



/ 6 



JL*. 



£ p - 2 



A, 

- y- 



0- SK')*' 



i 4 



r-f 



D 



X. 

-r 



Yti, z 



jr- ( 



/ rz' 






•Liepmann, H. W. , ^nd Pickett, A.E., "Introduction to Aerodynamics 
of a Compressible Fluid", John Wiley & Sons, 1947, p. 26. 
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whore the only variable for any one free stream is :r. 

The coefficient of drag, Cjj, is defined for a three dimensional 



bodyj 




I 



0 





z 



where d is the maximum diameter of the body. 



0 - 



J_ . . , >. 
2 /^V. lV 





c 



v 





w 

2. TT j\ Cp (Js\ 



This integral is numerically evaluated with C 



P 



being known from 



the above eauotion 
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DISCUSSION OF R2SULTS 

The Sauer Method of Characteristics was applied to the nose 
section of the Corporal 2 for free streara Mach numbers of 2 and 4. 



L - 96”, Tho result for a free streara Mach number of 3 was obtained 
from a solution for the 3ame nose section by H, K. Forster, For the 
case of M - 4 two entirely independent calculations were made to 
check the consistency of results. In one of these solutions the 
nose was assumed to be conical for a distance of six inches from 
the vertex and the constant velocities obtained from the Taylor- 
‘‘accoll solution for a cone were assumed to occur over this portion 
of the nose. In the second of the computations the nose was assumed 
to be conical for a distance of nine inches from the vertex. It 
can be seen from the plot of pressure coefficients obtained in Fig. 
(6) that the results were consistent. 

The pressure coefficients obtained for M - 2 and M * 4 are 
tabulated in Tables 17, V and VI and are plotted in Figs. (8) and 
(9). For comparison, values of Cp, obtained by the Karnan Moore 
Method as applied by V?. Z. Chien* to Ogive Nose Sections, are 
plotted. Fig3. (12) and (13). This method makes use of a linear- 
izing nrocesa and a close check of results would not be expected. 

It is noted however that a closer check is obtained at the higher 
Hach numbers. This is consistent with theory. 

Values of the drag coefficients obtained are! At K = 2, Cp- 
0.09168, at M = 3, CL, = .066, and at K = 4, 0 = 0.061. These are 
compared to the values obtained from the Karman Moore Method in 



The shape of the missile was given by the parabola 




Fig. (14). 



"■Wave Drag of a Projectile Nose at a Supersonic Velocity by the 
Karnan-Moore Method, ORDCIT Report Ho. 4-24, Anril 19, 1946. 
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All of the points computed in the Mach nets are indicated in 
Figs. (5), (6), end (7), end the velocities and shock wave angles 
are tabulated in Tables II and III. (In these tables the velocity 
is expressed in terms of the critical velocity, C", M* * ZQlg.gitff) 

The calcul £ tion of points and of velocities at these points 
theoretically could be carried out to any desired decree of accuracy 
defending on the scale used in the graphical solution. Actually 
the scale used is limited in the hodograph plane by the sice of 
drafting equipment available. In this instance a twenty-four inch 
rule calibrated in l/50 inch increments (5 units equal 1 inch) was 
used. In the M - 4 solution 50 units correspond to M* = 1, and in 
the <‘1 = 2 solution 65 units correspond to M* = 1. These scales 
wore chosen so thAt the full length of the rule would be utilised 
(and thus maximum accuracy) in the hodograph plane. In the physical 
plane 5/3 units corresnond to 1 inch (i.e. a 1:3 scale). 

Following is an outline of the steps in the solution of A = 4 
for the Corporal E nose section: 

1. Taylor-Maccoll solution* 

A. A value for U/c corresponding to M - 4 was chosen from the 
graph. 

B. A tabular form was set up and indicated calculations were 
made using increments of (ft equal to 0.5 degrees. 

C. The angle of the 3hock wave va3 solved for in the Te.ylor- 
Maccoll region as indicated. 

3. The next step was to solve for the exact free stream Mach 
number. 

•Taylor, G. I. end Maccoll, J. W. , "Air Pressure on a Cone Moving at 
High Speeds", Proc. of the Roy. Soc., A, Vol. 139, Fages 276-311. 1933, 
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S. Operations vero carried oat as indicated in Table VII 




rt / 
c ~ / 



J..Hi 7 



M* was converted to the scale in use 



2. The nose section was drawn to the scale chosen. 

3. A shock polar was constructed for Mach number 3.9712 as indicated 
in Fig. (5). 

4. Coordinates of the hodograph plane were drawn above the physical 
plane in such a position as not to interfere with the Mach net. 

5. A table of Mach angles and sine squared of the Mach angle was 
made for anticipated valocities. (See Table VIII). 



6. From Taylor-Maccoll 1 s solution velocities were drafted in the 



(Point 1 was taken a distance from the nose corresponding to 
6 inches in one case and 9 inches in the other). 

7. Calculations of the Mach net by the Sauer Method was started. 

It was advisable to calculate the slone of the nose at the inter- 
section by a J line from the geometry of the body rather than 
to approximate it graphically. 

8. When a point on the shock wave had to be found the shock polar 
was placed over the hodograph plane and was removed after the 
point was located. 

9. Z was calculated. (see Table IX). 

* 

10. 0^ was calculated, (see Table IV), and a plot of C^ vs. X was made. 

11. -hren increments of X were taken and the corresponding values of 
r and Cp were tabulated. 



* A. 




hodograph •''lane and in the physical plane the f) lines were 
drawn using average values of Mach angle between the points. 
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12, A plot rOp vs. r, see Fig. (10), vms made end tho ’'positive" 

area under this curve was determined. Thio area was connected 
to the proner scele. (In this case the conversion factor 
was (('$') * O'i )• 



13 The value for Cp = - ~~ x Area was calculated. 

It vas found that the most efficient method of constructing 
the Kach net was to work in pairs. One person drafted while the 
other performed the calculations.. From three to ten points per 
hour were completed after proficiency in the mechanical processes 
was obtained. The average working speed was about five points per 



hour. 



MACH NET FOR M = 3.9712 (Solution I) 
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MA.CH NET FOR M = 2.9712 (Solution II) 
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MACH NET FOR M = 2.0039 
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In Column 3, 

50 units correspond. 



to M* = 1 

r 



TABLE I 

3AMPL3 TABLE SHOVING DETAILED 
calculation 3 pop. *iagh hw 



p 


P 




cT 


* 


sin 2 * 


V 


69 


j69 

]69 - 68 
x 69 - 61 


102.05 

101.87 

102.05 
102.05 


2 

9 

9 

9 

9 


i 

15 

22 

39 

15 


2 

17 

17 


i 

9 

4 


.0870 

.0862 


16.35 

16.58 

17.22 

16.35 




l 70 


102.46 


8 


60 








15.87 




1 70 - 69 


102.25 


9 


2 


16 


57 


.0849 


16.11 




t70 - 62 


102.26 


9 


20 


16 


56 


.0849 


16.79 


70 




102.46 


8 


50 








15.87 




1?1 


103.12 


8 


15 








14.84 




l71 - 70 


102.79 


8 


32 


16 


38 


.0819 


15.35 




X 71 - 63 


103.22 


8 


60 


16 


24 


.0797 


15.96 


21 




103.16 


8 


S5 








15.10 




1?2 


103.88 


7 


58 








14.46 




l72 - 71 


103.57 


8 


11 


16 


04 


.0766 


14.78 


72 


* 


103.75 


7 


55 








14.44 




173 


101.54 


8 


44 








15.48 




l?o - 65 


101.37 


9 


14 


17 


27 


.0899 


16.10 


73 




101.54 


8 


44 








15.48 


74 


165 - 73 


101.37 


9 


7 








16.10 




175 


101.48 


8 


55 








15.82 




]7b - 74 


101.42 


9 


2 


17 


27 


.0899 


15.96 




i75 - 66 


101.40 


9 


27 


17 


27 


.0899 


16.32 


75 




101.50 


9 


6 








16.00 




176 


101.72 


8 


52 








15.80 




x 76 - 75 


101.61 


8 


58 


17 


20 


.0887 


15.90 




]76 - 67 


101.69 


9 


8 


17 


16 


.0881 


16.23 


76 


- 


101.72 


8 


52 








15.80 




177 


101.42 


9 


0 








16.04 




x 77 - 76 


101.57 


8 


56 


17 


21 


.0889 


15.92 




2?7 - 68 


101.56 


9 


15 


17 


21 


.0889 


16.43 


77 




101.73 


8 


50 








15.90 




j78 


102.00 


8 


48 








15.61 




,78 - 77 


101.87 


8 


49 


17 


10 


.0872 


15.76 




:78 - 69 


102.03 


9 


02 


17 


05 


.0864 


15.98 


78 


X 


102.00 


8 


48 








15.61 


Point 


73 is a point on 


the shock wave. 


Point 


72 is 


a point on the missile. 





All other points are free stream points. P denotes the column for estimated 
points. The individual estimated points are designated by ), g( )» etc. 

) - ( ) denotes average values between the indicated points. 
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TABLE I (Qont ' 


'd.) 












_ JP . 


r ...... 


A4- 


d n 




da 




ft 












( 


r 


f 




O 


.2 


TT 




22.40 




















22.65 


3.52 






.224 


26 


31 


-7 


47 




20.87 






.485 




26 


43 


-7 


25 


69 


22.40 




















21.83 




















22.11 


3.96 






.244 


25 


59 


-7 


55 




20.31 




6.90 


.483 




26 


16 


-7 


36 


70 


21.83 


< 


















20.90 




















21.36 


6.62 






.427 


25 


10 


-a 


6 




19.45 




6.82 


.445 




25 


14 


-7 


34 


71 


.. 20.90 




















19.96 




















20.43 


6.90 






.381 


24 


13 


-7 


53 


72 


19.96 




















30.34 




















27.18 




14.0 


.745 




26 


41 


8 


13 


73 


30.34 


















74 


27.18 




















26.82 




















26.95 


2.38 






.127 


26 


29 


-8 


25 




?5.27 




6.92 


.40 




26 


54 


8 


0 


75 


26.82 




















26.40 




















26.61 


2.83 






.15 


26 


18 


-8 


22 




24.89 




6.98 


.404 




26 


24 


-8 


8 


76 


26.40 




















26.00 




















26.20 


2.65 






.143 


26 


17 


-8 


25 




24.45 




6.90 


.412 




26 


36 


-e 


6 


77 


26.00 



















05.50 

25.75 3.75 .200 25 59 -8 21 

23.95 7.00 .403 26 7-83 

78 25.50 # 



TABL* II (a) 
3X0* V^OOI?H33 

In Column 2 'i = .5.9712 

50 units corresponds 3 :T I 

to > a = 1 



Aoint * O 

£ 1 



1 




98.60 


17 


05 


. 1 




98.65 


15 


44 


3 




98.81 


14 


38 


4 




99.05 


13 


17 


5 




9S.40 


12 


07 


6 




99.62 


11 


30 


7 




99.93 


15 


34 


8 




99.92 


13 


53 


9 




99.91 


15 


02 


10 




99.95 


12 


38 


11 




99.98 


13 


37 


12 




99.98 


14 


37 


13 




100.17 


11 


29 


14 




100.02 


12 


15 


15 




99.96 


13 


00 


16 




100.08 


14 


02 


17 




100.38 


10 


41 


18 




100.18 


11 


25 


19 




100.10 


12 


05 


20 




100.22 


13 


00 


21 




100.53 


13 


25 


22 




100.00 


11 


05 


23 




100.48 


10 


15 


24 




100.24 


10 


40 


25 




100.60 


10 


25 


26 




100.45 


10 


53 


27 




100.36 


11 


30 


28 




100.45 


12 


18 


2S 




100.74 


12 


40 


30 




101.17 


12 


50 


31 




100.82 


9 


48 


32 




100.61 


10 


30 


33 




100.56 


11 


0 


34 




100.55 


11 


40 


35 




100.79 


11 


53 


36 




101.11 


12 


Cl 


37 




101.25 


12 


17 


38 




101.30 


9 


10 


39 




101 . 06 


9 


29 


40 




100.79 


10 


0 


41 




100.65 


10 


34 


4° 




100.69 


11 


04 


43 




ICO. 88 


11 


13 


44 




101.17 


11 


17 


45 




101.37 


11 


36 



34 . 

FOE 



oint 




1 


4 






Si . 


i 


46 


101.79 


11 


29 


47 


101.07 


9 


50 


48 


100.94 


10 


10 


49 


100.94 


10 


45 


50 


101.03 


10 


55 


51 


101.41 


10 


50 


52 


101.51 


11 


10 


53 


102.00 


10 


55 


54 


102.50 


10 


49 


55 


101.30 


9 


10 


56 


101.18 


9 


43 


57 


101.08 


10 


07 


58 


101.27 


10 


15 


59 


101.51 


10 


14 


60 


101.65 


10 


15 


61 


102.05 


1C 


03 


62 


192.50 


3 


50 


63 


103.33 


9 


25 


64 


101.30 


9 


03 


65 


101.20 


9 


30 


66 


101.32 


3 


32 


67 


101.57 


9 


25 


68 


101.70 


c 


30 


69 


102.05 


c 

KJ 


15 


70 


102.46 


8 


50 


71 


103.16 


8 


25 


72 


103.75 


7 


58 


73 


101.54 


8 


44 


74 


101.37 


9 


07 


75 


101.50 


9 


05 


76 


101.72 


8 


52 


77 


101.73 


3 


50 


78 


102.00 


8 


48 


79 


ID 1 ’. 48 


8 


15 


80 


103.10 


7 


42 


81 


103.75 


7 


24 


82 


104.60 


6 


28 


83 


101.71 


8 


40 


84 


101.86 


8 


33 


85 


102.00 


r 


35 


86 


102.30 


8 


25 


87 


102.68 


8 


0 


88 


103.24 


7 


25 


89 ■ 


103.78 


6 


53 


90 


104.61 


r» 

o 


53 



91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 



oO| 



TABLES II (a) Cont'd. 



M* 


4 






2 


t 


105.24 


4 


50 


101.96 


8 


0 


101.83 


8 


20 


102.00 


8 


10 


102.15 


8 


15 


102.31 


8 


0 


102.70 


7 


40 


103.35 


7 


05 


103.85 


6 


25 


104.65 


5 


30 


105.30 


4 


38 


106.44 


2 


40 


102.13 


7 


48 


102.15 


7 


45 


102.43 


7 


48 


102.78 


7 


19 


103.41 


6 


35 


103.98 


6 


01 


104.75 


5 


03 


105.28 


4 


05 


106.33 


2 


21 


107.49 


0 


18 


102.45 


7 


20 


102.29 


7 


34 


102.44 


7 


25 


102.56 


7 


15 


102.93 


6 


53 


103.43 


6 


10 


103.98 


5 


30 


104.64 


4 


30 


105.25 


3 


35 


105.28 


1 


50 


107.25 


0 


12 


107.18 


0 


0 



*> 

‘ f 



In Column 
50 units corresponds 
to i " = 1 
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r 



TABLE II (b) 

FLOW VELOCITIES FOR M « 3.9712 
SH? II 



Point 


M* 




d 


Point 


M* 




Of 






£ 


i 






£ 


i 


1 


98.60 


17 


08 


46 


101.10 


11 


10 


2 


98.75 


15 


18 


47 


101.45 


11 


20 


3 


98.97 


13 


41 


48 


101.60 


11 


45 


4 


99.40 


12 


07 


49 


102.00 


11 


50 


5 


99.50 


16 


05 


50 


101.00 


S 


45 


6 


99.65 


13 


55 


51 


101.00 


10 


25 


7 


99.70 


15 


42 


52 


101.10 


10 


30 


8 


99.90 


18 


10 


53 


101.40 


10 


34 


9 


99.40 


14 


00 


54 


101.65 


10 


45 


10 


99.80 


15 


13 


55 


101.85 


10 


45 


11 


99.62 


11 


30 


56 


102.00 


11 


13 


1° 


100.10 


11 


25 


57 


101.25 


8 


50 


13 


99.55 


13 


00 


58 


101.12 


9 


18 


14 


99.85 


13 


50 


59 


101.10 


9 


45 


15 


ion. 00 


14 


47 


60 


101.25 


9 


50 


16 


99.86 


11 


28 


61 


101.50 


10 


00 


17 


100.15 


10 


30 


62 


101.60 


10 


15 


18 


99.85 


11 


45 


63 


101.90 


10 


20 


19 


100,00 


12 


15 


64 


102.00 


10 


40 


20 


100.00 


13 


15 


65 


102.60 


10 


24 


21 


100.15 


14 


12 


66 


101.30 


9 


15 


2? 


100.05 


10 


55 


67 


101.40 


9 


30 


23 


100.35 


10 


15 


68 


101.65 


S 


30 


24 


100.10 


11 


25 


69 


101.70 


9 


45 


25 


100.20 


11 


45 


70 


102.10 


9 


55 


26 


100.25 


12 


25 


71 


102.20 


10 


10 


27 


IOC;. 35 


13 


20 


72 


102.60 


9 


55 


28 


100.90 


13 


40 


73 


103.00 


9 


53 


29 


100.19 


11 


02 


74 


101.60 


8 


15 


30 


100.35 


10 


15 


75 


101.45 


8 


45 


31 


100.35 


10 


50 


76 


101.60 


8 


55 


32 


100.40 


11 


40 


77 


101.70 


8 


51 


33 


100.50 


12 


15 


78 


101.90 


9 


07 


34 


100.90 


12 


32 


79 


102.20 


8 


55 


35 


101.20 


13 


02 


80 


102.30 


9 


15 


36 


100.90 


9 


30 


81 


102.65 


9 


00 


37 


100.53 


10 


23 


82 


103.00 


9 


05 


38 


100.60 


10 


45 


83 


103.60 


8 


.30 


39 


100.60 


11 


35 


84 


101.80 


6 


25 


40 


101.00 


11 


45 


85 


101.95 


8 


20 


41 


101.40 


12 


0 


86 


102.00 


8 


50 


42 


101.50 


12 


25 


87 


102.25 


8 


45 


43 


101.72 


9 


56 


88 


102.40 


8 


50 


44 


100.90 


10 


12 


89 


102.70 


8 


40 


45 


100.90 


10 


55 


90 


103.05 


8 


30 



( 



91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

115 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 



TA3L2 II (b) Oont'd. 

r 



M * 


d 

& 


i 


103.60 


8 


10 


104.15 


7 


33 


102.10 


7 


35 


101.95 


8 


00 


102.00 


8 


10 


102.05 


8 


15 


102.30 


8 


05 


102.95 


7 


32 


103.95 


7 


00 


104.00 


6 


45 


104.80 


5 


50 


102.10 


7 


50 


102.20, 


7 


50 


102.50 


7 


42 


103.00 


7 


00 


103.55 


6 


15 


104.05 


6 


10 


104.85 


5 


15 


105.55 


4 


10 


102.30 


7 


00 


102.20 


7 


35 


102.40 


7 


30 


102.60 


7 


13 


103.10 


6 


40 


103.65 


6 


12 


104.10 


5 


45 


104.80 


4 


45 


105.50 


3 


50 


106.30 


2 


20 


102.35 


7 


10 


102.65 


7 


00 


103.20 


6 


25 


103.70 


5 


55 


104.10 


5 


25 


104.95 


4 


30 


105.45 


3 


38 


106.30 


2 - 


16 


107.30 


0 


12 


102.60 


6 


50 


102.48 


7 


00 


102.85 


6 


45 


103.30 


6 


20 


103.80 


5 


40 


104.35 


5 


10 


105.00 


4 


05 


125.65 


3 


20 


106.40 


2 


00 


107.40 


0 


05 


107.50 


0 


0 



Coluran 2, 
units corr 

?r= 1 


es-jonds 






TABLE 


III 
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FLO'i Vr&OCITISS 


P0R H = 2.0039 








Point 


M* 




cT 




Point 


M* 


< 


sT 






£ 




i 






SL 


• 


1 


94.90 


13 




46 


46 


99.38 


6 


13 


r. 

s. 


95.65 


11 




19 


47 


99.60 


o 


2C 


3 


96.38 


9 




47 


48 


99.05 


7 


00 


4 


97.28 


8 




22 


49 


99.70 


7 


55 


5 


98.01 


7 




25 


50 


98.70 


8 


30 


6 


99.44 


6 




55 


51 


98.65 


9 


10 


7 


98.94 


6 




00 


52 


98.90 


10 


20 


8 


99.19 


6 




05 


53 


99.00 


11 


00 


9 


96.05 


16 




03 


54 


99.69 


1? 


25 


10 


96.20 


12 




4C 


55 


99.58 


c 


56 


11 


95.90 


15 




49 


56 


99.95 


5 


50 


12 


96.65 


10 




35 


57 


99.40 


6 


35 


13 


96.20 


12 




35 


58 


99.10 


7 


20 


14 


95.85 


15 




30 


59 


99.20 


7 


A 5 


15 


97.40 


8 




55 


60 


99.10 


8 


15 


16 


97.00 


10 




10 


61 


99.10 


9 


5 


17 


96.30 


11 




58 


62 


99.15 


9 


30 


18 


96.50 


15 




00 


63 


99.70 


10 


30 


19 


97.30 


8 




00 


64 


100.35 


11 


43 


20 


97.G5 


8 




45 


65 


100.25 


5 


30 


21 


97.05 


10 




05 


66 


99.80 


6 


10 


22 


96.35 


12 




01 


67 


99.50 


6 


50 


23 


97.30 


14 




28 


68 


99.55 


7 


14 


24 


98.35 


7 




25 


69 


99.45 


7 


37 


26 


98.10 


8 




00 


70 


99.45 


8 


12 


26 


97.50 


9 




05 


71 


99.50 


8 


40 


27 


97.20 


10 




35 


72 


99.90 


9 


20 


28 


97.40 


12 




20 


73 


100.20 


10 


00 


29 


97.55 


14 




05 


74 


100.80 


11 


00 


3C 


98.75 


6 




50 


75 


100.80 


4 


35 


31 


98.65 


7 




12 


76 , 


100.43 


5 


12 


32 


98.15 


8 




00 


77 


100.10 


5 


40 


3o 


97.80 


9 




15 


78 


99.80 


6 


10 


34 


97.95 


10 




20 


79 


99.50 


7 


00 


35 


98.00 


11 




30 


80 


99.40 


7 


40 


36 


98.50 


13 




28 


81 


99.85 


8 


15 


37 


99.10 


6 




30 


82 


100.15 


8 


50 


38 


99,00 


6 




50 


83 


101.00 


9 


00 


39 


98.50 


7 




40 


84 


101.00 


10 


18 


40 


97.95 


8 




40 


85 


100.62 


4 


54 


41 


98.15 


9 




40 


86 


100.30 


5 


25 


42 


98.20 


10 




25 


87 


100.10 


5 


40 


43 


98,50 


12 




00 


88 


99.80 


6 


30 


44 


98.60 


13 




10 


89 


99.75 


7 


10 


45 


99.78 


5 




40 


90 


100.05 


7 


30 



91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 



r 



TABLE III (Cont'd.) 



M* 




of 


Point 


M* 




f 




SL 


1 






SL 


• 


100.30 


7 


50 


140 


102.00 


5 


10 


100.95 


8 


00 


141 


103.45 


4 


45 


100.30 


8 


55 


142 


104.45 


4 


30 


101.80 


9 


33 


143 


105.65 


4 


06 


100.50 


5 


05 


144 


101.20 


4 


30 


100.30 


5 


25 


145 


101.35 


4 


50 


100.20 


5 


45 


146 


101.80 


4 


55 


99.85 


6 


30 


147 


102.35 


5 


00 


100.10 


6 


50 


148 


103.75 


4 


30 


100.30 


7 


15 


149 


104.80 


4 


05 


101.20 


7 


10 


150 


105.90 


3 


40 


100.95 


8 


00 


151 


107.80 


r*s 


21 


101.55 


8 


20 


152 


101.24 


4 


08 


102.20 


8 


48 


153 


101.65 


3 


45 


100.75 


4 


30 


154 


101.44 


3 


56 


100.62 


4 


48 


155 


101.35 


4 


10 


100.50 


5 


05 


156 


101.40 


4 


25 


100.35 


5 


25 


157 


101.80 


4 


20 


100.30 


6 


00 


158 


102.35 


4 


20 


100.40 


6 


.20 


159 


103.70 


3 


45 


100.55 


6 


40 


160 


104.65 


3 


20 


101.35 


6 


35 


161 


105.70 


2 


50 


101.20 


7 


10 


162 


107.45 


1 


35 


101.65 


7 


30 


163 


108.20 


0 


46 


102.20 


7 


45 


164 


101.65 


4 


00 


103.05 


7 


57 


165 


101.75 


4 


15 


100.62 


4 


48 


166 


102.10 


4 


10 


100.45 


5 


05 


167 


102.60 


4 


05 


100.45 


5 


30 


168 


103.80 


4 


30 


100.50 


6 


00 


169 


104.70 


3 


00 


100.60 


6 


15 


170 


105.70 


2 


30 


101.40 


6 


00 . 


171 


107.40 


1 


12 


101.35 


6 


35 


172 


108.25 


0 


45 


101.65 


6 


50 


173 


109.00 


0 


00 


102.00 


7 


05 










102.80 


7 


05 










103.75 


7 


05 










101.25 


4 


05 










100.83 


4 


34 










100.75 


5 


10 










100.80 


6 


35 










101.40 


5 


50 










102.05 


6 


05 










103.60 


5 


50 










104.80 


5 


37 










100.04 


4 


20 










100.70 


4 


35 










100.80 


5 


00 










101.35 


5 


10 











?oii 

1 

7 

9 

12 

16 

21 

30 

37 

46 

54 

63 

72 

82 

91 

102 

112 

124 
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T Alii 2 IV 
DATA SHEET 
FOR CALCULATION 
OF C FOR M = 3.9712 
(CORRESPONDING TO FIG. (6). 



Coord 

& 


inate3 

A. 


50 H* 


M* 


°p 


‘' r0 p 


15 


4.36 


98.60 


1.9720 


.2109 


1.83905 


18.83 


5.48 


99.93 


1.9986 


.1621 


1.97662 


22.32 


6.45 


99.91 


1.9982 


.1634 


2.10786 


26.23 


7.49 


99.98 


1.9996 


.1604 


2.40279 


32.62 


9.15 


100.08 


2.0016 


.1596 


^ o 

C*m 2<oUOO 


37.96 


10.45 


100.53 


2.0106 


.1419 


2.96571 


43.19 


11.67 


101.17 


2.0234 


.1218 


2.84281 


48.74 


12.96 


101.25 


2.0250 


.1192. * 


3.08966 


55.48 


14.34 


101.79 


2.0358 


.1032 


2.95978 


63.05 


15.78 


102.50 


2.0500 


.0833 


2.62895 


75.20 


18.00 


103.33 


'’.0666 


.0619 


2.22840 


88.17 


19.96 


103.75 


2.0750 


.0519 


2.07185 


101.97 


21.75 


104.60 


2.0920 


.0330 


1.43550 


116.63 


23.12 


105.24 


2.1048 


.0201 


0.92 942 


136.40 


24.46 


106.44 


2.1288 


-.0016 


-0.78272 


157.52 


24.98 


107.49 


2.1498 


-.0178 


-0.88929 


182.10 


25.00 


107.18 


2.1436 


-.0133 


-0.66500 



i 
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TABLE V 

DA! TA SHEET FOR CALCULATION OF C FOR « = 3.9712 

P 

CORRESPONDING TO FIG. (5) 



oint 


X 


R 


50 M* 


M* 


C 

P 


1 


10 


3.03 


98.6 


1.972 


.2109 


5 


12.60 


3.80 


99.5 


1.980 


.1957 


7 


16.20 


4.75 


99.70 


1.994 


.1701 


10 


20.60 


6.00 


99.80 


1.996 


.1666 


15 


24.85 


7.10 


100.00 


2.000 


.1597 


21 


30.20 


8.45 


100.15 


2.003 


.1546 


28 


35.60 


9.75 


100.90 


2.018 


.1301 


35 


41.15 


11.10 


101.20 


2.020 


.1207 


42 


47.20 


12.50 


101.50 


2.030 


.1116 


49 


52.80 


13.80 


102.00 


2.040 


.0971 


56 


58.30 


14.95 


102.00 


2.040 


.0971 


65 


65.90 


16.40 


102.60 


2.052 


.0806 


73 


73.70 


17.70 


103.00 


2.060 


.0702 


83 


83.30 


19.30 


103.60 


2.072 


.0554 


92 


93.00 


20.65 


104.15 


2.083 


.0428 


101 


107.80 


22.30 


104.80 


2.096 


.0287 


109 


122.80 


23.70 


105.55 


2.112 


.0132 


119 


139.10 


24.60 


106.30 


2.126 


.0007 


128 


157.65 


24.85 


107.30 


. 2.146 


-.0151 


139 


180.80 


25.00 


107.50 


2.150 


-.0179 



rC 

P 



.6350 

.7436 

.8079 

.9396 

1.1333 

1.3063 

1.2684 

1.3397 

1.3950 

1.3400 

1.4516 

1.3218 

1.2425 

1.0692 

.8338 

.6400 

.3128 

.0172 

-0.3752 

-0.4475 



0 

9 

11 

14 

18 

23 

29 

36 

44 

54 

64 

74 

84 

94 

104 

116 

127 

136 

143 

151 

163 

173 



4 



TA3L3 VI 

DATA 3HB.TT FOR CALCULATION OF C p FOR M - 2.0039 
CORRESPONDING TO FIG. (9) 



X 


r 


65 M* 


K* 


°P 


*°T, 


10.00 


3.08 


94.57 


1.454922 


.28291 


0.871 


12.40 


3.80 


96.05 


1.477692 


.24400 


0.927 


16.00 


4.50 


95.90 


1.475384 


.24292 


1.121 


18.00 


5.35 


95.85 


1.474615 


.25051 


1.340 


22.80 


6.65 


96.50 


1.484615 


.23419 


1.557 


27.90 


7.95 


97.30 


1.496923 


.21442 


1.705 


31.50 


8.85 


97.55 


1.500769 


.20831 


1.844 


37.60 


10.30 


98.50 


1.515538 


.18513 


1.907 


40.50 


11.00 


98.60 


1.516923 


.18299 


2.013 


47.30 


12.60 


99.65 


1.533076 


.15822 


1.994 


53.80 


14.00 


100.35 


1.543846 


.14205 


1.989 


60.50 


15.40 


100.80 


1.550769 


.13178 


2.029 


67.00 


16.55 


101.00 


1.553846 


.12724 


2.106 


74.00 


17.80 


101.80 


1.566153 


.10934 


1.946 


80.80 


18.90 


102.20 


1.572307 


.10053 


1.900 


88.70 


20.10 


103.05 


1.685384 


.08205 


1.649 


96.40 


21.10 


103.75 


1.596153 


.06732 


1.420 


110.00 


22.60 


104.80 


1.612307 


.04528 


1.023 


123.60 


23.90 


105.65 


1.625384 


.02805 


0.670 


139.10 


24.70 


107.80 


1.658461 


-.01373 


-0.339 


153.20 


25.00 


108.20 


1.664615 


-.02122 


-0.531 


168.80 


25.00 


109.00 


1.676922 


-.03592 


-0.898 
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Table vii 

TAYLOH-MACCOLL 1 M = 2.0039 

VELOCITIES AT A MISSILE TIF FROM TAYLOR-! <AC COLL SOLUTION 



(t) 


-V 


u 


M* 




_ v; 


— 






* 


9 


r 








2 


i 


ic 


2 


i 


o l 


2 




17.08 


0 


.59700 


1.4549 


37 


12 


0 


< 


3 ~ 


17 5 


54 


17 


19.00 


.03852 


.59635 


1.4564 


37 


8 


.06459 


3 


42 


15 18 


52 


26 


21.00 


.07538 


.59435 


1.4600 


36 


56 


.12683 


7 


14 


lo 46 


50 


42 


23.00 


.10992 


.59111 


1.4651 


36 


44 


.18595 


10 


32 


12 28 


49 


12 


25.0 


.14283 


.58670 


1.4715 


36 


26 


.24344 


13 


41 


11 19 


47 


45 


26.5 


.16674 


.58265 


1.4768 


36 


15 


.28618 


15 


58 


10 32 


46 


47 


28.0 


.19017 


.57798 


1.4829 


35 


56 


.32902 


18 


13 


9 47 


45 


43 


29.5 


.21328 


.57210 


1.4892 


35 


39 


.37241 


20 


26 


9 4 


44 


43 


31.0 


.23626 


.56682 


1.4966 


35 


20 


.41682 


22 


38 


8 22 


43 


42 


32.0 


.25162 


.56257 


1.5019 


35 


8 


.44727 


24 


6 


7 54 


43 


02 


33.0 


.26713 


.55804 


1.5078 


34 


51 


.47869 


25 


35 


7 25 


42 


16 


34.0 


.28295 


.55325 


1.5144 


34 


34 


.51143 


27 


5 


6 55 


41 


29 


35.0 


.29938 


.54818 


1.5222 


34 


14 


.54613 


28 


38 


6 22 


40 


36 


36.0 


.31711 


.54280 


1.5319 


33 


50 


.58421 


30 


18 


5 42 


39 


32 


35.333, 


.30596 


.54639 


1.5260 


34 


06 


.55997 


29 


15 


6 05 


40 


10 



angle of shock wave 



0 = geometrical angle selected for Taylor-^ccoll solution 
(f) - arc tan 
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TABLS VIII 



VELOCITIES AT MISSILE TIP FROM TAYLOR-HACCOLL 

SOLUTION M = 3.9712 



<p 




u2 


H* 


IS 


_ yz 

7i^ 


O 


✓ 




e 


KL 

I 


13 


.00066 


.65448 


1.97253 


1% 04 


.031656 


1 


49 


16 11 


35 


15 


20 


.006117 


.651540 


1.97637 


18 67 


.096893 


5 


32 


14 28 


33 


25 


22 


.01656 


.64571 


1.98324 


18 44 


.16013 


9 


05 


12 55 


31 


39 


23.5 


.02764 


.63950 


1.99081 


18 27 


.20865 


11 


47 


11 43 


30 


10 


23.376 


















30 


21 


19 


.00427 


.65256 


1.97507 


18 59 


.06475 


3 


42 


15 18 


34 


17 


21 


.01072 


.64898 


1.97938 


18 51 


.12853 


7 


19 


13 41 


32 


32 


«*N . •* 

ciO 


.02371 


.64175 


1.98801 


10 36 


.19221 


10 


53 


12 07 


30 
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= geometrical angle selected for Taylor-Maccoll solution 



arc tan 
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